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1 Introduction 



Current knowledge of the Cabibbo-Kobayashi-Maskawa (CKM) matrix leads to the Stan- 
dard Model (SM) expectation that the mixing-induced CP violation phase in 5° decays 
proceeding via the b — )■ ccs transition is small and accurately predicted [l]. Therefore, 
new physics can be decisively revealed by its measurement. This phase denoted by 0s is 
given in the SM by — 2arg [V't^l^^/KsV^^], where the Vij are elements of the CKM matrix. 
Motivated by a prediction in Ref. [2], the LHCb collaboration made the first observation of 
5° J/ip fo{980), /o(980) vr+vr" fs], which was subsequently confirmed by others 
This mode is a CP-odd eigenstate and its use obviates the need to perform an angular 
analysis in order to determine 0^ (6], as is required in the J/ip(j) final state (7||8]. In this 
Letter we measure 0^ using the final state J/ipTi^ii^ over a large range of tt+tt" masses, 
775-1550 MeVQ which has been shown to be an almost pure CP-odd eigenstate We 
designate events in this region as /odd- This phase is the same as that measured in J/^p4' 



decays, ignoring contributions from suppressed processes 10 . 

The decay time evolutions for initial 5° and decaying into a CP-odd eigenstate, 
/_, assuming only one CKM phase, are RT] 



.,_) X ( ^Arst/2 g-Ar,t/2 
r (^5° ^ /_j =Afe'^'' < — — (l + cos0,) + (l-cos0s)±sin0,sin(Amst) 

(1) 

where AFg = Fl — Fh is the decay width difference between light and heavy mass eigen- 
states, Tg = (Fl + Fh)/2 is the average decay width, Am^ = mn — w^l is the mass 
difference, and A/" is a time-independent normalization factor. The plus sign in front of 
the sin 0s term applies to an initial and the minus sign to an initial B^ meson. The 
time evolution of the untagged rate is then 

r (5° ^ /-) +r ^ /-) =^^e-^^' |e^r^'*/2(l + cos0.) + e-^^^*/2(l-cos0.)| . (2) 

Note that there is information in the shape of the lifetime distribution that correlates AF^ 
and 4>s- In this analysis we will use samples of both flavour tagged and untagged decays. 
Both Eqs. [T]and[2]are invariant under the change 0^ — )■ tt — 0^ when AF^ — )■ — AF^, which 



gives an inherent ambiguity. Recently this ambiguity has been resolved 12 , so only the 
allowed solution with AF, > will be considered. 



2 Data sample and selection requirements 

The data sample consists of Ifb"^ of integrated luminosity collected with the LHCb 
detector |13j at 7 TeV centre-of-mass energy in pp collisions at the LHC. The detector is 
a single- arm forward spectrometer covering the pseudorapidity range 2 < rj < 5, designed 

^We work in units where c — h= 1. 
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for the study of particles containing 6 or c quarks. Components include a high-precision 
tracking system consisting of a silicon-strip vertex detector surrounding the pp interaction 
region, a large-area silicon-strip detector located upstream of a dipole magnet with a 
bending power of about 4 Tm, and three stations of silicon-strip detectors and straw drift- 
tubes placed downstream. The combined tracking system has a momentum resolution 
6p/p that varies from 0.4% at 5 GeV to 0.6% at 100 GeV, and an impact parameter (IP) 
resolution of 20 |^m for tracks with high transverse momentum (pt)- Charged hadrons 
are identified using two ring-imaging Cherenkov (RICH) detectors. Photon, electron and 
hadron candidates are identified by a calorimeter system consisting of scintillating-pad and 
pre-shower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons 
are identified by a muon system composed of alternating layers of iron and multiwire 
proportional chambers. The trigger consists of a hardware stage, based on information 
from the calorimeter and muon systems, followed by a software stage which applies a full 
event reconstruction. 

Events were triggered by detecting two muons with an invariant mass within 120 MeV 
of the nominal J/ip mass fl4||. To be considered a J/ip candidate, particles of opposite 
charge are required to have greater than 500 MeV, be identified as muons, and form a 
vertex with fit P^r number of degrees of freedom less than 16. Only candidates with 
a dimuon invariant mass between —48 MeV and +43 MeV of the J/ip mass peak are 
selected. For further analysis the four-momenta of the dimuons are constrained to yield 
the J/ijj mass. 

For this analysis we use a Boosted Decision Tree (BDT) [Ts] to set the J/'ipTc~^n~ 
selection requirements. We first implement a preselection that preserves a large fraction 
of the signal events, including the requirements that the pions have p^ > 250 MeV and 
be identified by the RICH. candidate decay tracks must form a common vertex that 
is detached from the primary vertex. The angle between the combined momentum vector 
of the decay products and the vector formed from the positions of the primary and the 
-Bg decay vertices (pointing angle) is required to be consistent with zero. If more than one 
primary vertex is found the one corresponding to the smallest IP significance of the 
candidate is chosen. 

The variables used in the BDT are the muon identification quality, the probability that 
the TT^ come from the primary vertex (implemented in terms of the IP x^)? the p^ of each 
pion, the vertex x^, the pointing angle and the B^ flight distance from production 
to decay vertex. For various calibrations we also analyze samples of B^ — > J/ipK*^, 
K*^ — 7r+fC~, and its charge-conjugate. The same selections are used as for J/ipn^iT' 
except for particle identification. 

The BDT is trained with 5^ — > J/'?/^/o(980) Monte Carlo events generated using 



Pythia 16 and the LHCb detector simulation based on Geant4 |17|. The following two 



data samples are used to study the background. The first contains J/ipn^T^^ and J/iI)t:~t:~ 
events with m{J/ilJ7T^7T^) within ±50 MeV of the B^ mass, called the like-sign sample. 
The second consists of events in the B^ sideband having m( J/t/'tt+tt") between 200 and 
250 MeV above the B^ mass peak. In both cases we require 775 < m(7r7r) < 1550 MeV. 
Separate samples are used to train and test the BDT. Training samples consist of 
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74,230 signal and 31,508 background events, while the testing samples contain 74,100 
signal and 21,100 background events. Figure [T] shows the signal and background BDT 
distributions of the training and test samples. The training and test samples are in 
excellent agreement. We select — )■ J/ipTi^T[~ candidates with BDT > to maximize 
signal significance for further analysis. 




BDT 



Figure 1: Distributions of the BDT variable for both training and test samples of J/^/'tttt 
signal and background events. The signal samples are from simulation and the background 
samples derived from data. 

The J/iI)t:^t:~ mass distribution is shown in Fig. [2]for the /odd region. In the 5° signal 
region, defined as ±20 MeV around the 5° mass peak, there are 7421±105 signal events, 
1717±38 combinatorial background events, and 66±9 rj' background events, corresponding 
to an 81% signal purity. The 7r"'"7r~ mass distribution is shown in Fig. [3j The most 
prominent feature is the /o(980), containing 52% of the events within ±90 MeV of 980 
MeV, called the /o region. The rest of the /odd region is denoted as /q. 



3 Resonance structure in the J/V^tt+tt final state 

The resonance structure in 5^ — )■ J/ipTx^-K^ decays has been studied using a modified 
Dalitz plot analysis including the decay angular distribution of the J /if) meson [9]. A 
fit is performed to the decay distributions of several tt+tt" resonant states described by 
interfering decay amplitudes. The largest component is the /o(980) that is described by 
a Flatte function [Ts]. The data are best described by adding Breit-Wigner amplitudes 
for the /o(1370) and /2(1270) resonances and a non-resonant amplitude. The components 
and fractions of the best fit are given in Table [TJ 

The final state is dominated by CP-odd S-wave over the entire /odd region. We also 
have a small D-wave component associated with the /2(1270) resonance. Its zero helicity 
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Figure 2: Mass distribution of the selected J/ipTi^-K^ combinations in the /odd region. The 
blue solid curve shows the result of a fit with a double Gaussian signal (red solid curve) and 
several background components: combinatorial background (brown dotted line), back- 
ground from B~ — )• J/tpK' and J/'ipTT~ (green short-dashed line), J/ipTT^Ti~ (purple 
dot-dashed), 5° — )• J/iprj' and 5° — )• J/'ipcp when — )■ 7r"'"7r~7r° (black dot-long-dashed), 
and — >■ J/ipK'n^ (light-blue long-dashed). 
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Figure 3: Mass distribution of selected tt+tt^ combinations shown as the (solid black) 
histogram for events in the B^ signal region. The (dashed red) line shows the background 
determined by fitting the J/tpn^n' mass in bins of tt+tt^ mass. The arrows designate the 
hmits of the /odd region. 
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Table 1: Resonance fractions in S° — )■ J/ipTx^n over the full mass range [9]. The final- 
state helicity of the D-wave is denoted by A. Only statistical uncertainties are quoted. 



Resonance 


Normalized fraction (%) 


/o(980) 


69.7 ±2.3 


/o(1370) 


21.2 ±2.7 


non-resonant tt+tt^ 


8.4 ± 1.5 


/2(1270), A = 


0.49 ±0.16 


/2(1270), |A| = 1 


0.21 ±0.65 



(A = 0) part is also pure CP-odd and corresponds to (0.49±0.16^q;o8)% of the total rate.j 
The |A| = 1 part, which is of mixed CP, corresponds to (0.21 ± Q.QbtQ^^)% of the total. 
Performing a separate fit, we find that a possible p contribution is smaller than 1.5% at 
95% confidence level (CL). Summing the /2(1270) |A| = 1 and p rates, we find that the 
CP-odd fraction is larger than 0.977 at 95% CL. Thus the entire mass range can be used 
to study CP violation in this almost pure CP-odd final state. 



4 Flavour tagging 

Knowledge of the initial 5° flavour is necessary in order to use Eq. [l] This is realized 
by tagging the flavour of the other h hadron in the event, exploiting information from 
four sources: the charges of muons, electrons, kaons with significant IP, and inclusively 
reconstructed secondary vertices. The decisions of the four tagging algorithms are indi- 
vidually calibrated using — t- J/ipK^ decays and combined using a neural network as 
described in Ref. 19 . The tagging performance is characterized by ^tag-D^, where et&g is 



the efficiency and D the dilution, defined as D = (1 — 2uj), where u is the probability of 
an incorrect tagging decision. 

We use both the information of the tag decision and of the predicted per-event mistag 
probability. The calibration procedure assumes a hnear dependence between the predicted 
mistag probability rji for each event and the actual mistag probability ui given by Ui = 
Po+Pi'iVi ~ iv)): where po andpi are calibration parameters and (77) the average estimated 
mistag probability as determined from the J/ip calibration sample. The values are po = 
0.392±0.002±0.009,pi = 1.035±0.021±0.012, and (r^) = 0.391. Systematic uncertainties 
are evaluated by using J/ipK^ separately from J/ipK" , performing the calibration with 

— 7- J/ipK*^ and B^ — )■ D*+/i^z7^ plus charge-conjugate channels, and viewing the 
dependence on different data taking periods. We find etag = (32.9 ± 0.6)% providing us 
with 2445 tagged signal events. The dilution is measured as D = 0.272 ± 0.004 ± 0.015, 
leading to £tag^^ = (2.43 ± 0.08 ± 0.26)%. 



^In this Letter whenever two uncertainties are given, the first is statistical and the second systematic. 
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5 Decay time resolution 



The 5° decay time is defined here as t = md ■ p/\p\'^ , where m is the reconstructed 
invariant mass, p the momentum and d the vector from the primary to the secondary 
vertex. The time resolution for signal increases by about 20% for decay times from 
to 10 ps, according to both the simulation and the estimate of the resolution from the 
reconstruction. To take this dependence into account, we use a double- Gaussian resolution 
function with widths proportional to the event-by-event estimated resolution. 



1 



=1 



(t-i-jHT 
2{Si<Tt)2 



(3) 



where t is the true time, at the estimated time resolution, fit is the bias on the time 
measurement, + /J = 1 are the fractions of each Gaussian, and Si and 5*2 are scale 
factors. 

To determine the parameters of T we use events containing a J/ip , found using a 
dimuon trigger without track impact parameter requirements, plus two opposite-sign 
charged tracks with similar selection criteria as for J/ipTc'^ir' events including that the 
J/i/jiT~^7r~ mass be within ±20 MeV of the mass. Figure |4] shows the decay time dis- 
tribution for this J/ipTT'^'K' prompt data sample for the /o region; the /o data are very 
similar. The data are fitted with the time dependence given by 



P™*W = (l-/i-/2)T(t;a,) + 



n 7-2 



T{t - i- at) 



(4) 



where /i and /2 are long-lived background fractions with lifetimes ti and T2, respectively. 
The resulting parameter values of the function T are given in Table [2j 



Table 2: Parameters of the decay time resolution function determined from fits to 
J/t/'tt+tt" prompt data samples. 



Parameter 


fo region 


/o region 


fit (fs) 


-3.32(12) 


-2.91(7) 


Si 


1.362(4) 


1.329(2) 


S2 


12.969(3) 


9.108(3) 




0.0193(7) 


0.0226(5) 



Figure [s] shows the at distributions used in Eq. |3] for J/ipn'^TT' events in the /odd 
region after background subtraction, and for like-sign background. Taking into account 
the calibration parameters of Table |2| the average effective decay time resolution for the 
signal is 40.2 fs and 39.3 fs for the /o and /o regions, respectively. The average of the two 
samples is 39.8 fs. 
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Figure 4: Decay time distribution of prompt J/ipn'^Ti candidates in the /o region. The 
dashed (red) hne shows the long-hved component, and the sohd curve the totaL 
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signal events after background subtraction, and for like-sign background. 



6 Decay time acceptance 

The decay time acceptance function is written as 



A(t;a,n,to) = C i^^^/°f,i n , (5) 
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where C is a normalization constant. The other parameters are determined by fitting the 
hfetime distribution of -B" — )■ J/'ipK*^ events, where K*'^ — )■ K~ti^ . Figure [6]^a) shows the 
J/ip K*^ mass when the K^tx^ invariant mass is within ±300 MeV of 892 MeV. There are 
155,743±434 signal events. The sideband-subtracted decay time distribution is shown in 
Fig. |6]^b) together with a hfetime fit taking into account the acceptance and resolution. 

1.82 ± 0.06, to = 0.105 ± 0.006 ps and a lifetime 



n 



This fit yields a = 2.11 ± 0.04 ps" 
of 1.516±0.008ps, in good agreement with the PDG average of 1.519±0.007ps 14 
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Figure 6: (a) Mass distribution of — > J/ipK*^ candidates. The dashed (red) line shows 
the background, and the solid (blue) curve the total, (b) Decay time distribution, where 
the small background has been subtracted using the mass sidebands. The (blue) curve 
shows the lifetime fit. 



We check our lifetime acceptance by comparing with a CDF measurement of the 
Bl J/ijfo effective lifetime of t"^ = 1.70l|J:i? ± 0.03 ps |5j obtained from a single 
exponential fitj^ A fit of the /o sample (see Fig. [7| yields r*^^ = 1.71 ± 0.03 ps, while we 
find T^^ = 1.67 ± 0.03 ps in the /o sample. These two values are consistent with each 
other, within the quoted statistical errors, and with the CDF result. 



7 Likelihood function definition 

To determine (ps an extended likelihood function is maximized using candidates in the B^ 
signal region 

£(0,) = e-(^-+^^^-) H P{m„ U, au, g,, rji) , (6) 

1=1 

where the signal yield, Ngig, and background yield, A^'bkg, are fixed from the fit of the 
J/ipTi^n' mass distribution in the /odd region (see Fig. [2|. A^obs is the number of 5° can- 
didates, rrii the reconstructed mass, U the reconstructed decay time, and au the estimated 
decay time uncertainty. The flavour tag, gj, takes values of +1, —1 or 0, respectively, if the 

■^This corresponds to the lifetime of the CP-odd eigenstate if (ps is zero (see Eq. pi. 
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Decay time (ps) 

Figure 7: Decay time distribution of — )■ J/ip /o candidates fitted with a single exponen- 
tial function multiplied by the acceptance and convolved with the resolution. The dashed 
line is signal and the shaded area background. 



signal meson is tagged as 5°, 5°, or untagged, and rji is the estimated mistag probability. 
Backgrounds are caused largely by mis-reconstructed 6-hadron decays, so it is necessary 
to include a long-lived background probability density function (PDF). The likelihood 
function includes distinct contributions from the signal and the background. For tagged 
events we have 

P{mi, ti, ati, qu r]i) = A^sig£:tagi'^^(mi)Pf ^(t^, g^, r]i\ati)P^f{ati) 

+A^bkgeta|P^'^^K)Pf ^(t,|a,jP:fs(a,,) , (7) 

where e^^^ refers to the fiavour tagging efficiency of the background. The signal mass 
PDF, P^s(m), is a double Gaussian function, while the background mass PDF, P^^{m), 
is proportional to e~"™ together with a very small contribution from P° — )■ J/iprj', A^^/, 
that is fixed in the (pg fit to 66 events obtained from the fit shown in Fig. [2j 

The PDF used to describe the signal decay rate, P^^^, depends on the tagging results q 
and 7]. It is modelled by a PDF of the true time t, R{t, q, if), convolved with the decay time 
resolution and multiplied by the decay time acceptance function found for P° — )• J/ip K*^ 
events. From Eq. [TJ it can be expressed as 

Q.1 v) e"'""* |cosh — + cos 0^ sinh — — q[l — 2uj{ri)] sin^^ sin(Amst) j> , (8) 

where uj{t]) is the calibrated mistag probability. Thus the PDF of reconstructed time is 
Pf^it, g, r]\at) = P(t, g, r]) ® T{t - t; at) ■ A{t; a, n, to) . (9) 
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For untagged events we use 



+iVbkg(l - e\ll)P:^^{m,)P^''\U\au)P^^^{au) . (10) 
The PDF describing the long-hved background decay rate is 

® T{t - i; at) ■ A{t; a'^'^^ n^'^^ tf^) , (11) 

where rj''^^, T2^^ and /g^'^^ parameterize the underlying double exponential function. The 
same functional form is used to describe the background decay time acceptance as for sig- 
nal (Eq. [5]) with different parameters that are determined by fitting the like-sign J/t/'tt^tt^ 
events in an interval ±200 MeV around the 5° mass. The P^'^((Tfj) and P^^^{ati) func- 
tions are shown in Fig. [5j The parameters that are fixed in the likelihood fit are listed in 
Table E 



Pt'-'itlat) 



1 - /2 

bkg 



t f 
^bkg ^ J2 



bke 



bkg 



bkg 
^2 



Table 3: Parameters used in the functions for the invariant mass and decay time describing 
the signal and background. These parameters are fixed to their central values in the fit 
for (/).,. 



Function 



Parameters 



P^s(m) 
P^kg(m) 



A^sig 



= 7421, A^bkg = 1717 ± 38, N^. = 66 ± 9 
mo= 5368.2(1) MeV, ^^"=8.1(1) MeV, a^=18.0(2) MeV, f^= 0.196(2) 
a = (-5.35 ± 1.15) X 10"^ MeV'^ 
f g ^ 0.65(5) ps, s = 2.0(8) ps, ^ = 0.06(2) 



^bkg ^ 3^22 



T(t — i; at) see Table [2 



(10) ps-\ n^'^g = 3.31(14), = 



ps. 



8 Results 

The likelihood of Eq. |6] is multiplied by Gaussian constraints on several of the model 
parameters. These are the LHCb measured value of Am^ = 17.63 ± 0.11 ± 0.02 ps~^ (20] , 
the tagging parameters po and pi, the decay time acceptance parameters to, a, and n, and 
both Ts = 0.657 ± 0.009 ± 0.008 ps"^ and AP, = 0.123 ± 0.029 ± 0.011 ps'^ given by the 
J/^jJ(j) analysis (T). The fit has been validated with full Monte Carlo simulations. 

Figure [s] shows the difference of log-likelihood value, Aln(£), compared to the one 
at the point with the best fit, as a function of (pg. At each value, the likelihood 
function is maximized with respect to all other parameters. The best fit value is 
0s = — 0.019!l^oj74]'^o;oo3 rad. (The systematic uncertainty will be discussed subsequently.) 
Values for 0^ in the /o and /o regions are — 0.26 ±0.23 rad and 0.29 ± 0.28 rad, respectively, 
consistent within the uncertainties. The decay time distribution is shown in Fig. |9j 
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Figure 9: Decay time distribution of 5° — )■ Jjip foAA candidates. The solid line shows the 
result of the fit, the dashed line shows the signal, and the shaded region the background. 



The presence of a sin 0s contribution in Eq. [T] can, in principle, be viewed by plot- 
ting the asymmetry [N{Bfj - iV(5°)] / [N{Bfj + N{B^)\ of the background-subtracted 



tagged yields as a function of decay time modulo 27r/Ams, as shown in Fig. 10 The 
asymmetry is consistent with the value of (ps determined from the full fit and does not 
show any significant structure. 

The data have also been analyzed allowing for the possibility of direct CP violation. 
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Figure 10: CP asymmetry as a function of decay time modulo 27t/ Arris- The curve shows 
the expectation for 0^ = —0.019 rad. 



In this case Eq. [8] must be replaced with 

_r f f , ArJ 2|A| , , AT J 

R{t, q,r]) oc e <^ cosh — h r-pj cos 0^ smh 

I 2 1 + |A|^ 



2 

q[l-2uj{r])] 



1 + IAI 



[2|A| sin0sSin(Am^t) - (1 - |AH cos(AmJ)] | . (12) 



The fit gives |A| = 0.89 ±0.13, consistent with no direct CP violation (|A| = 1). The value 
of (ps changes only by —0.002 rad, and the uncertainty stays the same. 

The systematic uncertainties are small compared to the statistical one. No additional 
uncertainty is introduced by the acceptance parameters. Arris, ^s, ATs or flavour tagging, 
since Gaussian constraints are applied in the fit. The uncertainties associated with the 
fixed parameters are evaluated by changing them by ±1 standard deviation from their 
nominal values and determining the change in the fitted value of (ps- These are listed 
in Table |4} The uncertainty due to a change in the signal time acceptance function is 
evaluated by multiplying A(t] a, n, to) with a factor (l+/3t), and redoing the B'^ — .J/ipK*^ 
fit with the 5° lifetime fixed to the PDG value. The resulting value of /3 = (1±3±3) x 10"^ 
is then varied by ±4.4 x 10~^ to estimate the uncertainty in (ps- An additional uncertainty 
is included due to a possible CP-even component. This has been limited to 2.3% of the 
total /odd rate at 95% CL, and contributes an uncertainty to (ps as determined by repeating 
the fit with an additional multiplicative dilution of 0.954. The asymmetry between 5° 
and production is believed to be small, and similar to the asymmetry between B^ and 



B^ production which has been measured by LHCb to be about 1% 21 . The effect of 
neglecting this production asymmetry is the same as making a relative 1% change in the 
tagging efficiencies, up for B^ and down for B^, which has a negligible effect on (ps- 
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Table 4: Summary of systematic uncertainties on (pg- Quantities fixed in the fit that are 
not included here give negligible uncertainties. The total uncertainty is found by adding 
in quadrature all the positive and negative contributions separately. 



v^uaniiLy 


A- AD 


+ Change — 


Change 






in (ps (rad) in (pg (rad) 


P 


A A ^ ^ 


0008 


-0.0007 




u.u^u 


-0.0006 


0.0014 


(ps) 


0.8 


-0.0014 


0.0014 


fbkg 


0.02 


-0.0006 


0.0012 


^bkg 


38 


0.0009 


-0.0001 
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0.0006 


0.0001 


mo (MeV) 


0.12 


0.0012 


-0.0004 


a"^ (MeV) 


0.1 


-0.0002 


0.0008 


a 


1.1 X 10-^ 


0.0003 


0.0003 


T function 


5% 


0.0005 


0.0005 


CP-even 


multiply dilution by 0.954 


-0.0008 




Direct CP 


free in fit 


-0.0020 




Total systematic uncertainty on (pg 


+U.004 
-0.003 





9 Conclusions 

Using 1 fb~^ of data collected with the LHCb detector, — )■ J/ip7i~^7i~ decays are 
selected and used to measure the CP violating phase (pg. The signal events have an 
effective decay time resolution of 39.8 fs. The flavour tagging is based on properties of 
the decay of the other b hadron in the event and has an efficiency times dilution-squared 
of 2.4%. We perform a fit of the time dependent rates with the lifetime and the 
difference in widths of the heavy and light eigenstates used as input. We measure a 
value of (ps = — 0.019lo;}74lo;Q03 rad. This result subsumes our previous measurement 
obtained with 0.41 fb~^ of data |6]. Combining this result with our previous result from 
B^s — > J/tlj(p decays [t] by performing a joint fit to the data gives a combined LHCb 
value of (pg = -1-0.06 ± 0.12 ± 0.06 rad. Our result is consistent with the SM prediction of 
— 0.0363;;'^q;qq}5 rad fll. In addition, we find no evidence for direct CP violation. 
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